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Abstract 

We present a configuration interaction method optimized for Fock-Darwin states of two- 
dimensional quantum dots with an axially symmetric, parabolic confinement potential subject 
to a perpendicular magnetic field. The optimization explicitly accounts for geometrical and dy- 
namical symmetries of the Fock-Darwin single-particle states and for many-particle symmetries 
associated with the center-of-mass motion and with the total spin. This results in a basis set of 
reduced size and improved accuracy. The numerical results compare well with the quantum Monte 
Carlo and stochastic variational methods. The method is illustrated by the evolution of a strongly 
correlated few-electron droplet in a magnetic field in the regime of the fractional quantum Hall 
effect. 

PACS numbers: 71.15.-m (Methods of electronic structure calculations), 73. 21. La (Quantum dots), 73.23.Hk 
(Coulomb blockade, single electron tunnelling) 
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INTRODUCTION 



In the configuration interaction (CI) method the Hamiltonian of an interacting system 
is calculated in the basis of many-electron configurations and diagonalized exactly. Unlike 
the correlated basis of Jastrow or Hylleraas functions, the CI basis functions do not au- 
tomatically include correlations among pairs of electrons, so their repulsive interaction is 
not minimized. Hence the size and choice of the basis influences the accuracy of results. 
Exact diagonalization (ED) has been an important tool used to invest iga te the electronic 
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251 ] . Not only does it provide benchmark results for the 
ground state energy and wave function, but also it gives access to excited states and hence 
permits an interpretation of a number of spectroscopic techniques such as far-infrared 



. and Raman 
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In contrast, the quantum Monte Carlo (QMC) 



mean-field Hartree-Fock \3l 



,0,0, HQ 



361 ] . and density- functional 



are restricted to the ground-state proper- 



theory (DFT) calculations 
ties. 

The objective of this work is to present a CI method optimized for the specific problem 
of interacting electrons in a parabolic dot in a magnetic field. The optimization explicitly 
accounts for geometrical and dynamical symmetries of single-particle states, and for many- 
particle symmetries associated with the center-of-mass (CM) motion and total spin. This 
results in a reduced basis size and improved accuracy, which in turn allows for a reliable 
computation of dot properties in the strongly correlated regime. Previous work on total spin 
and CM-resolved calculations used either the group-theoretical approach limited to a small 
number of electrons 0, Q], or the lowest Landau level approximation n]]. Here we build 
on these results to construct a reliable and versatile computational tool capable of testing 
many of the interesting predictions involving electronic correlations in quantum dots. The 
paper is organized as follows. In Sec. II we focus on the physical properties of our system, 
discuss the Hamiltonian, its symmetries and quantum numbers, and its scaling properties. 
In Sec. Ill we present a total spin- and CM-resolved CI technique. We compare our results 
with QMC and stochastic-variational method (SVM) and with previous results describing 
the evolution of the electronic droplet with the magnetic field. 
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HAMILTON! AN AND SYMMETRIES 



We consider a two-dimensional quantum dot with a parabolic confinement potential of 
strength uj (henceforth h — 1) in a perpendicular magnetic field B = (0, 0, B). The Hamil- 
tonian for iV particles in real-space representation [with r = (x, y), p = (p x ,p y ), and the 
vector potential A(r) = B/2(—y,x)] is: 

N / 1 1 \ 1 N p 2 
H = Y(—(i>j + eAfr,)) 2 + -m*uj 2 r 2 ) + - V ' r r, (!) 

with the Zeeman term omitted for simplicity. Here m*, e are the electronic effective mass 
and charge, respectively, uo c = eB/m* is the cyclotron energy, and e is the dielectric constant 
of the host semiconductor. We can rewrite the one-particle contribution of equation to 
obtain (pf/2m* + m*LU 2 ,r 2 /2 + (u c /2)lj\ where uJh = \JuJq + uj 2 /A is the hybrid frequency 
and lj is the angular momentum operator. Introducing the effective length l^ 2 = m*uj^ 
and scaling all lengths in the units of lh and energies in effective Rydbergs Ry = m*e 4 /2e 2 
results in a dimensionless Hamiltonian, 



■ (2) 



Here the strength of Coulomb interactions is proportional to the effective length (Ih/ao), with 
a = e/m*e 2 being the effective Bohr radius and Iq 2 = m*u c being the magnetic length. 
Therefore, the larger the effective length or the smaller the hybrid frequency, the larger the 
contribution from electron-electron interactions. 

The Hamiltonian ((TJ exhibits another scaling property. To see it we rewrite it as a 
function of u , u c , and N and separate the contribution of the angular momentum. This 
concept is similar to the modification of the Hamiltonian of a (natural) atom in Ref. 

The total Hamiltonian is therefore a sum of a Hamiltonian of a dot in zero magnetic field 



with the confinement frequency y lu 2 + uj 2 /4 } and the total angular momentum contribution 
{u c /2)L. As a result we can map the spectra and wave functions of a dot in the presence 
of the magnetic field to those at zero field. Note that Eq. Q implies that the energy levels 
of subspaces with different angular momenta are shifted against each other due to the term 
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(u c /2)L, whereas within a specific subspace of the angular momentum the energy levels are 
fixed relative to each other. Therefore, ground states in high magnetic fields are excited 
states in zero field and a suitable choice of a finite basis is not trivial. 

The choice of the basis is facilitated by the properties of the one- and the many-particle 
Hamiltonian. The properties of the one-electron Hamiltonian are easily understood by a 
transformation from the momentum and position representation (p, r) into the language of 
harmonic-oscillator raising and lowering operators 

H = u+(a\.a+ + 1/2) + w_(ala_ + 1/2), (4) 



represe 

m 



with the energy spectrum e nma = u) + {n + l/2)+oj_{m+l/2) (n, m — 0, 1, oo), and the two 
oscillator frequencies u± = Uh±uj c /2. The harmonic-oscillator spectrum is equivalent to the 
spectrum derived independently by Fock Q| and Darwin j^J (the Fock-Darwin spectrum, 
FD) applying different methods. The corresponding states are denoted by \nma), where 
a =|, I denotes the electronic spin. The eigenvalues of the angular momentum operator 
I = n — m, the spin s 2 , and its z-component s z are good quantum numbers. In addition, 
there are dynamical symmetries associated with the form of the confining potential. These 
symmetries lead to degeneracies of the harmonic-oscillator spectrum at B — and also at 
finite fields whenever uj + /uo^ = k, with k = 1,2, .... 

For numerical calculations it is useful to resort to the language of second quantization. 
The Hamiltonian of N particles in second quantization with Fermion creation $ nma and 
annihilation c nma operators can be written as: 

H = J2 £ iAaCia + r Yl (hj\W\k,l)cl a c] a/ c ka/ ci a , (5) 
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where the composite index i = (n, m) (similarly j, k, I) denotes the harmonic-osci 
quantum numbers. The Coulomb matrix elements (\W\) are given explicitely in Ref. 
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The Hamiltonian (0) conserves total spin, total z-component of spin, S z = 1/2 J2ia ^cj^c^, 
and total angular momentum, L = J2nma( n ~ m )Cn ma c nma . Apart from the spin and the 
geometrical and dynamical symmetries discussed above, the many particle Hamiltonian (j3J) 
las also a hidden symmetry related to the separability of the CM and relative motion 
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. The CM symmetry can be expressed by the introduction of new operators describing 
collective CM excitations of the system: 



N nil in 



it = 



-7= E Vm + lcl (m+1) 

V -i ' nmcr 



CT C-nma ; 



ramcr 

^+ = 7/^ E V / "'^(n-l)m<r^^' 



inner 



A- = -7= E v^c^m-i)^™- (6) 



nma 



The operators and At (i+ and A_) satisfy the following commutation relation with 
the full interacting Hamiltonian if: [H, A\] = u)±A± and [H, A±] = — u±A±. They can be 
interpreted as creation (annihilation) operators for collective excitations of the system since 
they increase (decrease) the energy of the system by uj±. They allow for the construction of 
exact eigenstates of the interacting system. 

We can now explicitly define CM operators C which commute with the interacting Hamil- 
tonian: 

C+ = A\A + — — \J ( n ' + ^) n c\n'+l)m'a'Crim'a'c\ n _ l)ma C nma , 

n'm'a'nma 
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C_ = AM_ = — VK + 1 ) m €>(m'+l)a'Cn>m>a>c\ m _ 1)lna Cnma- (7) 

n'm'a'nma 

Because Slater determinants constructed from a FD basis set are eigenstates of L and S z , 
it is natural to resolve the total angular momentum L and the total spin in z-direction S z 
as good quantum numbers. On the other hand, the structure of the operators S 2 and C± 
in second quantization is similar to that of the Coulomb interaction, involving two creation 
and two annihilation operators. In what follows we construct ED methods which make use 
of these operators to reduce the size of the Hamiltonian matrix. 



SPIN-RESOLVED CI METHOD 



The form of the total spin operator in the language of creation and annihilation operators 

is: 

S 2 = N/2 + 4 2 -Es f T 4wT- (8) 

Its last term annihilates one particle in a state i and recreates it on the same orbital, but 
with flipped spin; the same is true for the orbital j. Thus, the pattern of orbital quantum 
numbers is conserved, which makes it possible to arrange the Slater determinants in blocks 
which enclose all states with the same pattern, and to rotate the single blocks into eigenstates 
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of S 2 . The size of the blocks can be obtained using combinatoric arguments. We start with 
a given electron number N, among which electrons are spin up, and spin down, and 
we distribute them on M orbitals (M < N). Because of Pauli's exclusion principle we get 

N 

— < max (N h N{) <M < N. (9) 

Thus we have N — M doubly occupied orbitals (singlets), M — N± unpaired spins up, and 
M — iVf unpaired spins down. The distribution of the singlets is determined by Pauli's 
principle. The number of different orbitals available for unpaired spins up and unpaired 
spins down is M — (N — M). The size of the block equals the number of all possible 
distributions of unpaired electrons spin up and down on M — (N — M) orbitals: 

(M-{N - M)\ (M-(N-M)-(M- N{)\ 



blocksize - 

(2M — N)l 



M-N i J \ M - iV T 



(M — JV;)!(M — iVf)!' ^ 
In Table H] we show the number of total spin eigenstates per block for a six-electron system as 
a function of M and S z = (iV| — iVj)/2. These blocks are fairly small (e.g., blocksize = 252 
for N = M — 10, S z = 0) and the diagonalization of the blocks has to be done only once 
for each value of M. 

In Table |H] we show how the size of the Hamiltonian is reduced by applying the scheme 
on a system with N = 5, L = 1, S z — 1/2 and 352954 Slater determinants and on a system 
with iV = 6, L = S z = and 326120 Slater determinants. Note that for even electron 
numbers N the largest subspace is that of S = 1, not S = 0. In the case of odd N the 
S — 1/2 subspace is larger than the others. 

We now apply the spin-resolved CI method to the calculation of the ground state energy 
of 4 through 7 electrons. The accuracy of the method is assessed by comparison of our 
calculations in zero magnetic field with available results from quantum Monte Carlo and 
stochastic-variational method. The calculations were performed using GaAs parameters 
(e = 12.4, m* = 0.067) and a confinement energy ujq = 3.32 meV. Table IIHI shows the 
ground state energies for a number of subspaces L, S. We see that the QMC and SVM 
energies appear to be slightly lower then the ED results, however, their order appears to 
change depending on the actual subspace and number of electrons. Altogether, the ED 
appears to give highly reliable results for ground and excited states which can be extended 
to finite magnetic fields. 
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TABLE I: Number of total spin eigenstates per block for a six-electron system as a function of M 
and S z . 



N 


S=l/2 


S=3/2 


S=5/2 


E 


5 


184769 


137705 


30480 


352954 



N 


5 = 


5 = 1 


5 = 2 


5 = 3 


E 


6 


92410 


152460 


70711 


10539 


326120 



TABLE II: Effective size of the Hamiltonian for systems with N = 5, L = 1, S z = 1/2 and 352954 
Slater determinants and with N = 6, L = S z = and 326120 Slater determinants. 



N 


(L,5) 


Eqmc [meV] 


E'svMlmeV] 


-E'ED [meV] 


rel. error of ED 


4 


(0,1) 


44.0439 


44.0261 


44.0763 


< 1.2 • 10~ 3 


4 


(-2,0) 


44.5301 


44.4945 


44.5498 


< 1.3 • 10~ 3 


4 


(0,0) 


44.8301 


44.8004 


44.8355 


< 8 • 10~ 4 


5 


(-1,1/2) 


65.6159 


65.5827 


65.6488 


< 10~ 3 


6 


(0,0) 


90.1166 


90.1391 


90.2383 


< 1.5 • 10" 3 


7 


(-2,1/2) 


118.9784 




119.1426 


< 1.2 • 10" 3 


7 


(0, 1/2) 


119.3045 




119.4069 


< 1.0 • 10~ 3 



TABLE III: Comparison of ground state energies for = 4, 5, 6, 7 of the subspaces L, S (uj = 3.32 
meV) at B = 0. QMC energies are from Ref. [3^, and SVM energies from Ref. j^. SVM results 
for seven electrons are not available. 
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SPIN- AND CM-RESOLVED CI METHOD 



The ED method with both the total spin and the CM resolved as good quantum num- 
bers is particularly applicable in the high magnetic field strongly correlated regime, where 



correlation effects dominate the reconstruction of the dot through magic states 



spin textures 



2l| , and skyrmions 10 1 . 



The construction starts with defining the number of electrons N and their total angular 
momentum L. Moreover, we restrict the single-particle basis set by choosing the number 
of Landau levels of interest. This is done by fixing the maximal quantum number n max . 
Since the total angular momentum is defined using the FD numbers n and m, the above 
choices impose a limitation on the number m, i.e., determine an m max such that J2f=i m j = 
J2jLi n j — L < n max N — L < m max . The Hilbert space is thus finite and we can generate all 
configurations and diagonalize them with respect to the CM and total spin operators. Since 
we restrict ourselves to n max Landau levels, only the CM operator C_ is of interest. From its 
definition (Eq. (JIJ)) we see that C_ increases/decreases the quantum number m, but does not 
change the pattern of numbers no. We can therefore organize our configurations in blocks 
with the same no pattern; these blocks are diagonalized separately. The general scheme is 
thus similar to that for the total spin, only here the CM diagonalization has to be done for 
each block as the matrix elements of (7_ in second quantization explicitly depend on m. We 
retain only the eigenstates with eigenvalue C_ = 0, because the eigenstates corresponding 
to eigenvalues 1, 2, 3, .. can be generated systematically by applying the CM motion creation 
operator it on eigenstates with eigenvalue 0. 

Our next aim is to resolve the total spin in the basis of stored CM eigenvectors. We have 
established that S 2 couples states with the same nm pattern. It is therefore convenient to 
introduce superblocks with the same n pattern. Thus, a superblock encloses all blocks with 
the same no pattern and S 2 can couple states only within a superblock. After calculating all 
eigenstates with zero CM motion in the previous step, all states with the same n pattern are 
collected and diagonalized with respect to the total spin operator S 2 . Note that one set of 
CM eigenvectors calculated for the system with lowest possible spin polarization (0 for even 
number of electrons, 1/2 for odd number of electrons) can be used to calculate all possible 
eigenvalues of S 2 . As a result we obtain a basis set consisting of eigenstates of the operators 
L, S z , S 2 , and C*_, which we use to generate the Hamiltonian matrix. Table UVl shows how 



S 



the Hilbert space is reduced for the (S z , L) = (0, —18) subspace of a four-electron system in 
the two-Landau- level approximation (n max = 1). The first column labels the superb locks. 
The second column contains the information which electrons are in the lowest (n = 0), and 
which are in the second Landau level (n = 1), and how many configurations are possible 
(third column). Each superblock may enclose several blocks which are labelled in the fourth 
column. The ncx-pattern and the corresponding number of configurations are in the fifth 
and sixth columns. After applying the CM operator and filtering out states with zero CM 
motion the number of states is reduced for each no block (seventh column), and to 1555 in 
total (see last row). In the next step all superblocks are diagonalized with respect to the 
total spin operator. The resulting numbers of configurations for S = 0, 1, and 2 are given in 
the last column. Thus, instead of having to solve an eigenvalue problem of the size 15456 in 
the subspace S z = 0, L = —25 we end up with three eigenvalue problems with matrix sizes 
528, 791, and 236 (se e last row). To illustrate the method and compare it with previous 
work we show the results for four electrons in the range of magnetic 

fields — 8 T for GaAs parameters (e = 12.4, m* = 0.067), confinement energy ujq = 3 meV, 
and one and two Landau levels. In Fig. ^ we show the energies of the ground and the two 
lowest-lying excited states of each subspace L, S with C_ = 0, without the Zeeman energy 
(g = 0). The spectrum in Fig. ^exhibits an energy band consisting of ground states of all 
S'-subspaces, is separated from the excited states by AE ~ 0.6 meV (marked by the double 
arrow). Since the energy needed to construct an excitation of the CM motion is of the same 
order of magnitude (c<j_ = 0.62 meV for 8 T and higher for smaller magnetic fields), this 
separation clearly shows that the lowest-lying excited states of the quantum dot for magnetic 
field between 4 and 8 T are all excitations of the relative motion. 

In Fig. El we show the evolution of the ground state energy with the magnetic field in 
the one- and two-Landau-level approximation with GaAs Zeeman energy (g=-0.44). In the 
one-Landau-level approximation we do not find the Hund's rule state expected at B = 
nor the v — 2 spin-singlet droplet with S — 0, L — —2. The first quantum Hall state is the 
spin-polarized state with S = 2, L = —6. Starting at B = 4 T we find the reconstruction 
of the v = 1 droplet. This correlation-driven effect is expressed in a number of ground 
state transitions with decreasing angular momentum (Fig. Efa), (b)). The ground states 



change both angular momentum and spin, in agreement with previous work 



□ □ 



In Fig. Efb) we show the evolution of the ground state in the two-Landau-level 
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Superblock 


n pattern 


Size 


Block 


n<7 pattern 


Size 


C_ = 


Spin resolved 


1 


(0,0,0,0) 


728 


a 


(o |,o |,o r,o T) 


728 


78 


5 = 0: 26 
5 = 1: 42 
5 = 2: 10 


2 


(0,0,0,1) 


3458 


a 
b 


(o 1,0 1,0 |,i t) 
(o |,i |,o |,o t) 


1729 
1729 


182 
182 


5 = 0: 126 
5=1: 182 
5 = 3: 56 


3 


(0,0,1,1) 


5880 


a 
b 

c 


(oi, oi, i|, it) 
(oi,i 1,0 T, i T) 
(l i.i |,o T,o T) 


910 
5880 
910 


91 
406 
91 


5 = 0: 196 
5=1: 301 
5 = 2: 91 


4 


(0,1,1,1) 


4270 


a 
b 


(o i, 1 1, 1 1, i T) 
(i I,i l,o T, i T) 


2135 
2135 


210 
210 


5 = 0: 145 
5=1: 210 
5 = 2: 65 


5 


(1,1,1,1) 


1120 


a 


(U,U,iT,iT) 


1120 


105 


5 = 0: 35 
5 = 1: 56 
5 = 2: 14 


S 




15456 






15456 


1555 


5 = 0: 528 
5=1: 791 
5 = 2: 236 



TABLE IV: Reduction of the Hilbert space for the (5 2 ,L) = (0,-25) subspace of a four-electron 
system in the two-Landau- level approximation using CM and total spin as good quantum numbers. 

approximation. We still miss the Hund's rule B = state, but recover the v = 2 state with 
L = — 6 and 5* = 0. 

The qualitative agreement between the lowest- (Fig. Efa)) and two-Landau- level 
(Fig. Efb)) approach is good for high magnetic fields, i.e., beyond the v = 1 quantum 
Hall droplet. The only difference is the (S, L) = (0, —8) phase which is present in Fig. Ufa) 
, and missing in Fig. Efb). The absence of the low-spin phase is due to the weakening of 
Coulomb interactions by inclusion of higher Landau levels. The second effect is that the 
fields corresponding to the ground-state transitions are lower in the one-Landau-level due 
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to the fact that in this approximation the correlation effects are underestimated. 

An advantage of the CI method over QMC and SVM methods is the ability to calculate 
excited states. In Fig. |3] we show the energy gap between the ground and the lowest-lying 
excited states as a function of the magnetic field and Zeeman energy for the N=4 electron 
system, calculated in the two-Landau-level approximation. The magnetic field changes from 
from to 10.5 T, where the filling factor v = 1/3-droplet (angular momentum —18) emerges. 
In Fig. Eta) we plot the energy gap between the ground state and the lowest excited state 
with the same S z (selection rule AS Z = 0, without restrictions for the quantum numbers 
L, S). This selection rule applies to charge excitations. Consequently, we cannot expect a 
systematic increase of the gap as a function of the Zeeman energy. All the structures in the 
diagram for increasing Zeeman energy are caused by transitions to higher polarized states. 
In contrast, we used the selection rule AS Z = ±1 for Fig.^b). This kind of excitation should 
be observed in microwave-radiation or in spin-flip-Raman experiments. The characteristic 
features in this case are energy gaps linearly growing with Zeeman energy. This behavior is 
especially pronounced in the regime of the v = 1 droplet around B = 3 T and in the regime 
of the v = 1/3 droplet around B = 10 T, where both ground and first excited state are 
stable. 



CONCLUSIONS 



We presented a configuration interaction method optimized for Fock-Darwin states of 
two-dimensional quantum dots with an axially symmetric parabolic confinement potential 
subject to perpendicular magnetic field. The optimization explicitly accounts for geometrical 
and dynamical symmetries of Fock-Darwin single-particle states and for many-particle sym- 
metries associated with the center-of-mass motion and with total spin. As a key element we 
introduced a block structure to cluster Slater determinants which are coupled by S 2 and/or 
C_. By diagonalizing these small blocks we managed to rotate the many-particle basis set 
into eigenstates of S 2 and/or C_, and used them to calculate the Hamiltonian matrix. This 
reduced the size of the basis set and increased the accuracy of exact diagonalization calcula- 
tions for quantum dots. The results compare well with available quantum Monte-Carlo and 
stochastic-variational results. Thus we show that the CI method gives reliable results for 
the ground and excited states and allows to study the effects of electron correlation and the 
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magnetic field over a broad range of system parameters. It is intended serve as a benchmark 
for other methods which can be extended to larger particle numbers. 
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FIG. 1: Energies of the ground state and the two lowest-lying excited states of four electrons 
with angular momenta L = — 25 to calculated in the two-Landau-level approximation. The 
confinement energy is u = 3 meV and the Zeeman energy is neglected. All these states are 
eigenstates of the center-of-mass operator with the quantum number C_ = 0; solid lines correspond 
to total spin for S = 0, dashed lines - to S = 1, and dotted lines to S = 2. 
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FIG. 2: Phase diagrams for a four-electron system with confinement energy 3 meV and GaAs 
Zeeman energy (g = —0.44) as a function of the magnetic field from to 8 T in one-Landau-level 
(a) and two-Landau- level approximation (b). 
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FIG. 3: Energy gap between the ground and the lowest-lying excited state with the same S z (a) 
and with AS Z = ±1 (b) for a four-electron system as a function of magnetic field and Zeeman 
energy. 
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